Statistical Hadronization of Supercooled Quark-Gluon Plasma 



(N 
O 

o 

>>'. 

S 



A. Keraner|] and J. Manninenjj] 

Department of Physical Sciences, 
P.O. Box 3000, FIN-90014 University of Oulu, Finland 

L. P. Csernai^j 

Section for Theoretical Physics, Department of Physics, 
University of Bergen, Allegt. 55, 5007 Bergen, Norway and 
KFKI Research Institute for Particle and Nuclear Physics, 
P.O.Box 49, 1525 Budapest, Hungary 

V. Magasg 

Center for Physics of Fundamental Interactions ( CFIF), Physics Department, 
Instituto Superior Tecnico, Av. Rovisco Pais, 1049-001 Lisbon, Portugal 

The fast simultaneous hadronization and chemical freeze out of supercooled quark-gluon plasma, 
created in relativistic heavy ion collisions, leads to the re-heating of the expanding matter and to the 
change in a collective flow profile. We use the assumption of statistical nature of the hadronization 
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with different quark-gluon plasma equations of state. 
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I. INTRODUCTION 

The hydrodynamical models have been used exten- 
sively to study the evolution of the hot, strongly inter- 
acting matter created in relativistic heavy ion collisions. 
These models apply on the space-time region, where the 
initial, hard processes have reached the stage where the 
local thermal equilibrium can be assumed, and the strong 
interactions between constituent particles are very fre- 
quent. 

On the other hand, the thermal statistical models ap- 
plied to describe the final hadron abundances have been 
very successful for nucleus-nucleus (AA) jl], ||, @, |), 
proton- nucleus (pA) ^ and even for elementary pp, pp 
and e + e~ [Q, ^ reactions. The latter ones, particularly, 
suggest the statistical nature of the hadronization pro- 
cess. 

In this work, we study the fast hadronization and 
chemical freeze-out (CFO) of locally thermalized quark 
gluon plasma (QGP) described by hydrodynamical evo- 
lution. This process is idealized by sudden hadronization 
over a three dimensional hypersurface. The matter cross- 
ing the hypersurface is controlled by conservation of en- 
ergy momentum- and relevant current conservation laws, 
and additionally, by assumption of apparent thermal and 
chemical equilibrium distribution of the resulting hadron 
gas (HG) . By apparent equilibrium distribution we mean 
that although there is clearly no room for kinetic equili- 
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bration in elementary collisions - or in fast, simultaneous 
phase transition and CFO in nuclear collisions - the final 
hadron spectra are dictated by maximum entropy bound 
by conservation of energy and charge densities. 

It is shown, that the first order quasistatic phase tran- 
sition is far too slow |)| |h]] to give the hadron chemical 
decoupling time of icFO ^5 10 fm/c determined exper- 
imentally, see [pd] and references therein. Thus, in or- 
der to avoid the entropy decrease in the non-equilibrium 
hadronization, QGP must be allowed to supercool, i.e. 
to develop mechanical instability before the phase tran- 
sition |pd| . Although the simultaneous phase transition 
and CFO are assumed to be a non-equilibrium process, 
we can still exploit the apparent equilibrium parametriza- 
tion - the consequence of statistical nature of hadroniza- 
tion. 

We will show that the shock-like hadronization of su- 
percooled QGP leads to the change in collective flow pro- 
file, and re-heating of the system. We study the system 
with various choices of the equation of state (EoS) on 
both sides of the FO hypersurface. 



II. THE CHEMICAL FREEZE OUT PROCESS 

The fast hadronization process, idealized to take place 
on the zero-volume hypersurface, leads to a discontinu- 
ity in energy-momentum- and charge conservation equa- 
tions. In general, this leads not only to change in density 
quantities in LRF, but to the change in flow velocity pro- 
file as well 0,0. 

Denoting post freeze out quantities by subscript FO, 
the energy-momentum (T) conservation at the surface 
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element da v leads to a set of four equations 



(1) 



In the isospin symmetric case, the electric charge conser- 
vation is trivial, and we are left with local baryon number 
and strangeness conservation equations: 



Ngda u = Ng FO da v 



N v s da v 



(2) 
(3) 



Here, we restrict our considerations to time-like FO 
hypersurfaces in order to avoid currents entering from 
post-FO to pre-FO side. The complications in the gen- 
eral case, where the space-like FO is included, are dis- 
cussed extensively in |l2|, For the time-like hyper- 
surface, there is always a proper Lorentz transformation 
for each point x such that solving the equations (0j2|^) 
can be carried out in a frame, where flow velocities are 
u{x) — 7(1, 0, 0, v) and u-fo(x) — 7fo(1j 0, 0, i>fo) on pre- 
and post-FO sides, respectively. After solving the FO 
equations, the actual flow velocity ufo = Tfo^wfo) 
is obtained by a simple spatial rotation. Provided with 
these tools, only two of the four equations (|l|) are inde- 
pendent. 

Usually, the strangeness in hydrodynamical simula- 
tions is assumed to be homogeneously distributed, so the 
condition n${x) = ns^o{x) — binds the parameters 
associated with strangeness conservation. 

The process described is generally non-adiabatic, but 
the entropy constraint 



(4) 



must be satisfied in every point on the FO hypersurface. 

A. The Statistical Post Freeze Out Hadronization 

Statistical models were applied for high energy colli- 
sions |f|, 15, |l6) from the beginning of this field. In the 
last decade a significant development of these models and 
the extension of the area of their applicability took place. 
The main reason for this is a surprising success of the sta- 
tistical approach in reproducing new experimental data 
on hadron multiplicities in nuclear (AA) 0,0, [| ||, ||] 
and elementary (e + e~, pp, pp) collisions H pfl. One of 
the important results of the analysis of hadron yield sys- 
tematics at high energies (SPS and higher) done within 
the statistical models is the approximate independence 
of the temperature parameter (T = 160 -7- 180 MeV) of 
the system size and collision energy. Another important 
feature found in statistical analyses is the common freeze- 
out condition, energy per hadron E/N ~ 1 GeV, found 
for all systems from GSI 2 GeV A NiNi to SPS 160 GeV 
A PbPb collisions |l7j . These results can be attributed 
to the statistical character of the hadronization process. 

Calculations within statistical models are straightfor- 
ward when the mean number of particles of interest is 



large and consequently it is enough to fulfill the conser- 
vation laws in the average sense, i.e. the grand canonical 
(GC) description can be used. In GC ensemble, the mean 
particle multiplicities are just proportional to the volume 
V of the system. Thus, the particle densities and the ra- 
tios of the multiplicities of two different particle species, 
which are usually used for the comparison with the data, 
are volume independent. 

This simple volume dependence is however not valid 
any more for a small system in which the mean particle 
multiplicity is low, like in pp collisions. In this case the 
material conservation laws should be imposed exactly on 
each charge configuration of the system, i.e. canonical 
ensemble (CE) description should be used. This condi- 
tion introduces a significant correlation between particles 
who carry conserved charges (see |l^, |l9| , and references 
therein). 

In PbPb collisions at the SPS collider and AuAu reac- 
tions at RHIC we produce large and hot systems, where 
the corrections are negligible and GC description is satis- 
factory for most of the produced particles. Nevertheless, 
the number of strange particles created is rather low, 
and the net strangeness is exactly zero. Therefore the 
CE description would be preferable for the strange par- 
ticles (we are not going to include particles with charm 
or more heavy quarks). However, the basic assumption 
of the hydrodynamics is the local thermal and chemical 
equilibrium, whereas the canonical conservation laws can 
only be realized globally. This prevents us using CE in 
post FO, which is nothing but the extension of the hy- 
drodynamical part with different EoS. This drawback is 
not of a serious concern, because the conditions reached 
in high energy heavy ion systems are expected to fulfill 
the requirements for GC description JlQj] . 

In the relativistic GC ensemble for ideal hadron gas 
(HG), the properties of the matter in unit volume 
are parametrized by temperature T and fugacities \ 
(or chemical potentials fj,i) for conserved charges. For 
example, the pressure and energy density are P = 
P{T, Ab, As) and e = e(T, A#, As) if we omit the fugacity 
Xq associated with electric charge conservation. 

In case of sudden freeze out there might not be suf- 
ficient time to achieve chemical equilibration of number 
of strange- and anti-strange quarks. Thus, an over- or 
underpopulation of strange hadron species may persist 
after hadronization. This can, however, be treated in the 
GC approach as we will see. 

Now that the thermal ideal gas of hadrons provides us 
with tools to calculate post-FO quantities, we can col- 
lect the set of equations needed to describe the simulta- 
neous FO and hadronization taking place at a time-like 
hypersurface. For the sake of simplicity, we choose a hy- 
persurface of constant coordinate time, do v = (1, 0, 0, 0), 
which is, of course, connected with any time-like hyper- 
surface by a proper Lorentz transformation. Recalling 
the energy-momentum, baryon current and strangeness 
conservation equations yields 

(e + P) 7 2 -P = (e FO + Pfo)7fo - Pfo (5) 
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(e + P)w 7 2 



(e FO + -Pfo)"fo7fo ( 6 ) 

nB.FOlFO (7) 

™.S\fo7fo = 0. (8) 



Given the pre-FO quantities, the solution of this 
set of equations gives the post-FO parameters: 
wfo: 2fO) A^fo and ^s,fo, flow velocity of the fluid ele- 
ment, temperature and fugacities for conserved currents, 
respectively. These parameters describe completely the 
hadron spectra and other statistical quantities at the 
LRF of the fluid element. 

In order to take into account the possible over- or 
underpopulation of strange- and anti-strange quarks in 
post-FO side, a parametrization of the conservation of 
the number of s, s pairs must be introduced. In the GC 
formulation, a conservation law gives rise to a fugacity 
parameter. The one associated with the number of s, s 
pai rs, n StS , is usually called 75, first introduced in Ref. 
p0[ . The assumption of the survival of n StS over the FO 
hypersurface yields yet another equation to be included 
in the set above: 



5 7 = (n s , s ) F0 7FO- 



(9) 



III. NUMERICAL STUDIES 

For quantitative studies of the FO, we have chosen a 
framework of Bjorken model |2l| for the expanding QGP. 
This allows us to cover many FO scenarios with relative 
ease. Within the Bjorken model, the evolution of mat- 
ter in one spatial dimension z and charge densities rii is 
governed by equations 



de 



P 



T 

dP 
dy 
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0, 



(10) 

(11) 
(12) 



where r = \lt 2 — z 1 is the proper time and y is the rapid- 
ity of a given fluid element. For vanishing net charges, 
thermodynamical quantities are constant along constant 
proper time curves on (<, z) plane. Taking the FO to 
take place at constant coordinate time, the proper time 
of freezing out fluid element decreases with increasing 
spatial coordinate, so the temperature is an increasing 
function of z. The constant coordinate time choice is 
made in order to have a clear picture of the FO process 
in different thermal circumstances. In addition to the 
equations above, one needs an equation to describe the 
thermodynamical structure of the evolving matter, the 
equation of state. 



A. The Equation of State - Supercooled QGP 

Three different equations of state for supercooled 
quark-gluon plasma are considered. The QGP is assumed 
to be an ideal gas of three flavors (u,d,s) of quarks and 
their antiquarks. We will model the thermodynamics of 
the QCD vacuum in terms of the MIT bag model |2^| by 
introducing a phenomenological bag constant B. In this 
model, confinement can be thought of as a rule which al- 
lows quarks to inhabit only regions of a " false vacuum" , 
which has minimum energy density B. This is due to the 
fact that quarks can not be found in physical vacuum 
but only inside hadrons that are small bags of the false 
vacuum containing triplets or pairs of quarks. In QGP, 
all the space is endowed with this energy density B. 

Light quarks (u,u,d,d) can be considered as massless, 
but the mass of strange quark can not be neglected. For 
massless quarks and gluons, integral in partition function 
Zqgp can be evaluated analytically, but the contribution 
of strangeness must be calculated numerically. It is con- 
venient to split Zqqp into four parts: 

In Zq GP = In Zf + In Z b + In Z s + In Z vac , 

where lnZ/, InZf, and hiZ s are the contributions of 
massless quarks, gluons and massive strange quarks re- 
spectively. In Z vac is the contribution of the bag constant 
added to the energy-momentum tensor. Evaluation of 
analytical parts of In Zqgp leads to the expressions |2^| 



lnZ f 
\nZ b 

In Z vac 



BV 
- — > 



1 



1 



(13) 
(14) 
(15) 



where fi q is a light quark chemical potential, [i q = /x_b/3. 
The total strangeness is zero during the whole collision 
evolution. If one assumes that the strangeness produced 
in a collision is spread homogeneously in QGP, then ns 
must be zero in any given fluid element, so the contribu- 
tion of s and s quarks is 



In Z s 



6V 

7T 2 



dpp In ( 1 



(16) 



where m s is the mass of the strange quark, here chosen 
to vary between 150 and 250 MeV. 

If one puts the bag constant to zero (B=0), then In Z vac 
term vanishes and we end up with the normal quantum 
distribution for the ideal gas. In this case, the explicit 
form of the three flavor QGP EoS is[f|§ 

P(T,(j,i) = -e, e(T,m) = e SB . (17) 

In the case of finite bag constant, EoS reads: 

P{T,fx B ,B) = -e S B-B, e(T,fi B ) = e S B+B, (18) 
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The third equation of state in addition to ones with 
finite and zero MIT bag constants considered here was 
suggested in the Ref. p 4 B5 | (and first time applied for 
heavy ion collisions in p6| ) for the case of baryon free 
QGP (n q = 0) with two quark flavors: 

P{T)= l -e{T)-bT, e(T)=e SB . (19) 

We call this a spinodal EoS, inspired by the existence of 
a local minimum in pressure profile (|l9|). It is easy to 
check that e and P are connected by the thermodynam- 
ical equation 



£ = T dT 



(20) 



This EoS can be justified as resulting from the non- 
perturbative QCD effects: The lattice calculations (see 
for example ^?J) show that close to the critical tem- 
perature for hadronic matter, QGP phase transition en- 
ergy density has a sharp increase and soon saturates with 
an equilibrium Stefan-Boltzmann value e(T) = esB{T), 
while pressure increases much more slowly and reaches 
Stefan-Boltzmann limit only at very high T. Such a pres- 
sure suppression is presented by eqs. (fl9|). 

The spinodal EoS considered here will be of form (|l9| ) 
but with three quark flavors and finite net baryon den- 
sity: 



l 



e{T,m)-bT, e(T, tH )=e S B. (21) 



Comparing these expressions with the bag model EoS 
one can see that for the pressure the bag constant, £?, is 
replaced by T and fi dependent term bT, while for energy 
density the analogy with bag model disappears due to the 
relation (20). 

If the phase transition temperature T c is known, 
parameter b can be fixed from the Gibbs condition 
Phg(T c , fi) — Pqgp(T c , /i). This same condition fixes 
the constant B within the bag model. 



B. Freeze Out Illustrated 

The post-FO matter is described by the quantum ideal 
gas, composed of hadrons up to a mass of 2.5 GeV, 
listed in the year 2000 issue of Review of Particle Physics 
pgj . In order to satisfy the conservation equation (g) for 

strange quark pairs, additional fugacity 75 is introduced 

\s- 1 

for each strange particle i as Xi — > Xijg ■ Additionally, 
the mesons carrying ss pairs must be taken into account. 
In this work, the number of meson i is affected by the 
fugacity factor = jg Cs , where c s is the relative ss con- 
tent in the meson. We take c s = 0.5 for 77 mesons, and 
c s = 1 for 0, / (980), /i(1510), 0(1680) and </> 3 (1850). 

In figure |l| we show the pressures of the HG and QGP 
as functions of temperature. The parameters appearing 
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FIG. 1: Hadron gas pressure and quark-gluon plasma pres- 
sures with ideal gas, bag model and spinodal model equations 
of state with parameters fiB = 100 MeV and T c = 160 MeV. 



in different equations of state for QGP are fixed by set- 
ting the critical temperature to T c = 160 MeV. The pres- 
sure minimum in the spinodal EoS is labeled by T m . Of 
course, the critical temperature varies with baryon den- 
sity (or hb)i but we find this variation to be negligible 
within reasonable range of FO density. The strangeness 
saturation parameter 75 is set to one, corresponding to 
the full strangeness equilibration. We always let the QGP 
side be in full strangeness equilibrium, so the Gibbs con- 
dition for the adiabatic, isothermic phase transition com- 
pels 75 = 1. In the following, however, the phase transi- 
tion is neither isothermic nor adiabatic, so the change in 
fugacities is allowed and unavoidable. 

The figure || depicts the variation of some relevant 
variables along the coordinate time t = 10 fm/c in the 
Bjorken model. Going from the edge of the light cone 
(z = 10 fm) to the non- flowing center of the system 
(z = fm), we describe the changes of quantities in 
proper time of the given fluid element from r = to 
t = 10. This is, going from the origin to the right, we 
go from the later, cooler and more dilute stage of evo- 
lution to the earlier, hotter and more dense stage. We 
have chosen the initial values for energy- and baryon den- 
sity eo = 5 GeV fm -3 and (713)0 = 1 fm -3 , respectively, 
at the time r = 1 fm/c. For the ideal gas QGP EoS, 
the higher eo = 8 GeV is used in order to keep the fig- 
ure more descriptive. The model parameters B in the 
MIT bag EoS and b in the spinodal EoS are fixed to pro- 
duce critical temperature 160 MeV, while for the ideal 
gas EoS there is no critical temperature. The arbitrary 
value T c = 160 MeV is chosen to give similar scaling for 
temperatures in all cases. The ratio R of entropy currents 
in HG and QGP stays well above one for bag and spin- 
odal models all the way up to unrealistic FO conditions, 
so the substantial increase in entropy is expected in fast 
simultaneous FO and hadronization. For the reference, 
the R for the ideal gas QGP is, as expected, far below 
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FIG. 2: Various quantities on the FO surface. Panels from 
top to bottom: a) ideal gas EoS for QGP, b) MIT bag EoS, 
c) spinodal EoS. R is the ratio of entropy currents in HG 
and QGP, Th and Tq are the temperatures on the HG and 
QGP side, and T c is the critical temperature. The model 
parameters in the bag- and spinodal EoS are fixed to produce 
critical temperature 160 MeV. For the ideal gas EoS, there is 
no real T c , the chosen 160 MeV is arbitrary. The vertical line 
in the panel c marks the point, where the HG temperature 
crosses the critical value. 



FIG. 3: Various quantities on the FO surface. Panels from 
top to bottom: a) ideal gas EoS for QGP, b) MIT bag EoS, c) 
spinodal EoS. vh and vq are the HG and QGP flow velocities. 
7Si and 7s 2 are the strangeness saturation fugacities with s 
quark masses 150 MeV and 250 MeV, respectively. The model 
parameters in the bag- and spinodal EoS are fixed to produce 
critical temperature 160 MeV. For the ideal gas EoS, there is 
no real T c , the chosen 160 MeV is arbitrary. The vertical line 
in the panel c marks the point, where the HG temperature 
crosses the critical value. 
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one. The crossing of HG temperature Tjj with T c marks 
the endpoint of the physically allowed FO. For the bag 
model, there is no such crossing, but the FO stays un- 
physical with all reasonable values of QGP temperature 
Tq. We could choose lower initial value for the QGP en- 
ergy density, but the the degree of supercooling, Tq/T c , 
would be of order 0.5 at the HG critical point. For the 
spinodal model QGP, the case is different. At the point 
Th/T c = 1, the degree of QGP supercooling is about 
24%, and stays rather constant when going towards the 
cooler system. 

In figure || we illustrate the changes in flow veloci- 
ties and strangeness saturation parameter in the circum- 
stances equal to ones in figure |[ For both cases, bag- 
and spinodal QGP, the flow is decelerated at the FO to 
HG, indicated by the ratio of HG and QGP velocities, 
vh/vq. At the point Tjj = 160 MeV for the spinodal 
EoS, this ratio is 0.80, indicating the final HG flow ve- 
locity of 0.54c. The 75, resulting from the survival of 
ss pairs through the FO process, is calculated for two 
different values of s quark mass m s . jg 1 corresponds to 
rn s = 150 MeV and js 2 is f° r m s = 250 MeV. Comparing 
the figures, we find 7s to be very sensitive to the choice 
of QGP EoS and the m s . It is worth noting, that fixing 
7s to one or varying the m s gives only negligible change 
to other quantities in figures and ^. 

IV. CONCLUSIONS 

We have found that the realistic and accurate study 
of the freeze-out process is important and cannot be ne- 



glected. Our results show that the FO process is very 
sensitive on the properties of the EoS. Since our main 
goal is to identify the EoS from the data this is an obser- 
vation of basic importance. 

The correct treatment shows that FO is not even pos- 
sible from arbitrary kind of initial state, and the entropy 
constraint is a sensitive condition. 

From the point of strangeness, we can also conclude 
that strangeness production is very sensitive to the 
correct FO treatment and to the pre FO EoS. Thus, 
strangeness data not only provide a signal of QGP forma- 
tion, but with proper and realistic description of freeze 
out, strangeness provides the most sensitive signal indi- 
cating different properties of the pre FO EoS. 

In conclusion, this study shows that collective, contin- 
uum reaction models, like fluid dynamical (FD) models 
(one fluid FD, multi - fluid FD, chiral FD), can and must 
be supplemented with realistic Freeze Out treatment to 
evaluate measurable data. These calculations indicate 
that this is now possible cell by cell in FD models, al- 
though it requires more than average computational ca- 
pacity, and needs preferably high performance parallel 
computing. 
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